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As a step towards deriving universal amplitude ratios of the superconductive phase transition 
we calculate the vacuum energy density in the symmetric phase of 0(iV)-symmetric scalar QED in 
. D = 4 — e dimensions in an e expansion using the minimal subtraction scheme commonly denoted 

by MS. From the diverging parts of the diagrams, we obtain the renormalization constant of the 
vacuum Zv which also contains information on the critical exponent a of the specific heat. As a 
side result, we use an earlier two-loop calculation of the effective potential []lj to determine the 
«4 ■ renormalization constant of the scalar field Zj, up to two loops. 
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Ctf ! I. INTRODUCTION 
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One of the most intriguing problems in the physics of critical phenomena is a theoretical understanding of the 
superconductive phase transition within the renormalization group approach. A first discussion was given in 1974 by 
Halperin, Lubensky, and Ma on the basis of the Ginzburg-Landau or U(l) Abelian-Higgs model, in 4 — e dimensions, 
' generalizing a similar four-dimensional analysis of Coleman and Weinberg j3|. In a one- loop approximation, they 
did not find an infrared-stable fixed point, and in spite of much effort it is still unclear whether a higher-loop 
renormalization group analysis would be capable of explaining the existence of a critical point in 4 — e dimensions. 
O , Experimentally, this existence has been confirmed only recently with the advent of high-T c superconductors. In 
conventional superconductors, the Ginzburg criterion or the more relevant criterion for the size of phase fluctuations 
fNJ \ H > predicted a too small temperature interval for the critical regime to observe anything beyond mean-field behavior. 
Evidence had so far come only from Monte-Carlo simulations |6| and an analogy with smectic-nematic transitions in 
' liquid crystal 0. Only by artificially allowing for an unphysically large number of replica n of the complex field cj>, 
larger than 365, has it been possible to stabilize the renormalization flow in 4 — e dimensions. Historically, a dual 
disorder formulation of the Ginzburg-Landau model brought in 1982 the first theoretical proof for the existence of a 
tricritical point at a Ginzburg parameter w 0.77, a material parameter characterizing the ratio between magnetic 
and coherence length scales ||, |t| • This prediction was recently confirmed by extensive Monte-Carlo simulations [jl0| . 
The confusing situation in the Ginzburg-Landau model certainly requires further investigation in higher loop approx- 
i imations. So far, two loop renormalization group calculations in 4 — e dimensions have not yet produced satisfactory 
results |ll], |l^]. Analyses in d — 3 dimensions a la Parisi have also left many open questions flifl . An interesting 
observation was made by Nogueira p"4j ], that an anomalous momentum instability below T c may be responsible for 
the unusual resistance of the superconductive transition to theory. Hope for a better understanding has also been 
raised by a recent renormalization group study in d — 3 dimensions performed for the first time below T c where a 
fixed point has been found at the one-loop level fl5|| . 

Once an infrared-stable fixed point is located, it will determine critical exponents and amplitude ratios. The former 
can be extracted from the perturbative expansions of the renormalization constants corresponding to coupling, mass 
and wave function renormalization. The latter also requires the calculation of finite parts of certain quantities. For 
instance, the amplitude ratios for the specific heat can be obtained by computing both the divergent and finite parts 
' of the vacuum energy in the symmetric and in the symmetry-broken cases. 

Our work will provide analytic results for the divergent as well as finite parts of the vacuum energy up to three 
loops in the symmetric case. While the divergent part is the same in the symmetry-broken case, the determination of 
the finite part in that case is left for future work. 

We work with the so-called modified minimal subtraction scheme, denoted by MS Jll| . The singularities are collected 
in the dimensionlcss renormalization constant Z v . In the course of the calculations, we also recover information about 
the other renormalization constants of the theory. This may be viewed as a cross check of our work. 

Since the theory has so far only a fixed point for large n > 365, we shall keep an arbitrary number of replica in the 
theory, the physical case being n = 1. The n complex fields are coupled minimally to an Abelian gauge field which 
describes magnetism. The N — 2n real and imaginary parts of the fields are assumed to have an 0(-/V)-symmetric 
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quartic self-interaction. 



II. MODEL 



The Lagrangian density to be studied contains n = N/2 complex scalar fields 
potential Ab^ and reads, with a covariant gauge fixing, 
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coupled to the magnetic vector 
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where Db^ — <9 M — iesAsfi denotes the covariant derivative, -Fb m „ 



duAsfj, is the field strength, and a a 



gauge parameter. The bare character is indicated by the subscript 'B.' In principle, there are also ghost fields which, 
however, decouple in the symmetric phase and remain massless. Working in dimensional regularization they do not 
contribute to the energy because of Veltman's rule, 
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The coefficient 1/4 in front of the coupling constant gs is conventional. The Feynman diagrams associated with 
the vacuum energy of the Lagrangian have been generated iteratively in Ref. [[l7| . At some places it will be useful 
to compare our results with those of an earlier work , where we have derived the two-loop effective potential above 
and below T c . For such comparisons, a replacement qb — > 2gs/3 is required. 

A full extension of the work in Ref. |l| is highly nontrivial since the effective potential requires the calculation 
of Feynman diagrams with three different masses. For this reason we shall restrict ourselves in this paper to the 
symmetric phase T > T c , where the field expectations vanish and the system contains only two masses, which 
greatly simplifies the problem, in particular since one of the masses, the photon mass, is zero. As a consequence, most 
diagrams can be reduced to scalar integrals which can be computed exactly. The only exception is the watermelon — or 
basketball — diagram whose e expansion is, however, known to sufficiently high order in e |fis|| . 

As in Ref. @], we shall use throughout Landau gauge a — > 0, which enforces a transverse photon field. This has 
the advantage of being infrared stable 20 . 



III. RENORMALIZATION 



The renormalization constants of the model are defined by 
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where, in the last equation, we have taken into account the relation Z e = Z<p, which is a consequence of the Ward 
identity. Heuristically, this equality comes from the requirement that the covariant derivative Db^b should not only 
be invariant with respect to gauge transformations but also with respect to renormalization. Thus it must acquire the 
same normalization factor as the field itself, going over into sj r Z^D /J/ (p. An arbitrary mass scale /x in Eq. (|^) serves to 
define dimensionless coupling constants g and e. 

The above multiplicative renormalizations are not sufficient to extract all finite information from the theory. The 
vacuum energy requires a special treatment, as emphasized in a previous work of one of the authors (B.K.) 
Dimensionality requires the effective potential to have mass dimension D. To have a finite vacuum energy we must 
add to the Lagrangian a counterterm 
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The different renormalization constants may be expanded in powers of the fluctuation size H as follows: 
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where the subscript j stands for fields and coupling constants <f>, A,m, g,e. In minimal subtraction, each expansion 
coefficient has simple power series 

zf = 9 l - k e 2k (c). k e- 1 + c 2 ^ 1 + ... + c',^ 1 ) , (6) 

k=0 

except for Z^ , where the systematics is 
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Initially, one finds also pole terms of the form 1/e 2 x In, 1/e x In, and 1/e x In 2 , where In is short for \n(m 2 /ji 2 ) with 
\x being related to the mass scale \i via the Euler-Mascheroni constant je as ft 2 = Airfi 2 exp(— je)- These, however, 
turn out to cancel each other, which provides us with a nice consistency check of the renormalization procedure p2[ . 

IV. FEYNMAN RULES AND VACUUM DIAGRAMS 

The elements of the Feynman diagrams associated with the Lagrangian (Q) are 
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In Ref. |l71| , the vacuum diagrams of the theory have been generated recursively up to four loops. Table 1 shows all 
diagrams needed for the three-loop vacuum energy above T c , omitting those which contain massless separable loop 
integrals which vanish by Veltman's rule (0). 



RENORMALIZATION CONSTANTS Z ( v s) AND 



The determination of the two-loop effective potential below T c JjJ allows us to extract the following one- and two- 
loop contributions to the renormalization constants. Note that a factor (47r)~ 2i has been taken out in the definition 
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A. Results for the Feynman integrals 

In this section, we give the value of the diagrams listed in Table 1. Although the exact value of part of the three-loop 
diagrams is known, we only give the s expansion through order s , for the sake of brevity. The notation is as follows: 
is the integral for the case N = 2, omitting the weights of Table 1 and setting coupling constants and scalar mass 
equal to unity. The first index I is the loop order while the second index n counts through the diagrams within each 
loop order as listed in Table 1. 
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where IA± in Eq. ( |25| ) denotes the polylogarithm of order 4: "Li^ix) = J2k=i x /k P 



B. Divergent terms through three loops 

Denoting by E v&c the vacuum energy density in the symmetric phase, we find up to three loops the expansion 
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A few remarks are useful on the calculation of the pole terms of the diagrams. Taking into account the expansion 
of the renormalization constants in the form (JsJ) , each Feynman integral is expanded up to order h 3 . The expansion 
coefficients are determined to cancel the 1/e™ terms arising from the Feynman integrals. In particular we have: 

- To order h, the cancellation of the 1 /e pole leads directly to the known one-loop result (|l8|) as in Ref. (jj . 

- To order ti 2 , there are pole terms of the form 1/e 2 , 1/e and 1/ex In. The 1/ex In terms cancel if 
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which is fulfilled by the one-loop expressions ( |l4|) and (|l6) for Z^l and Z^\ respectively. After this, the 
cancellation of the ordinary poles lets us recover the result ( 2l]) obtained before in Ref. B . 

Finally, we need to cancel the poles at the h 3 level. Besides poles without logs, there are poles of the types 
1/e x In, 1/e 2 x In and 1/e x In 2 , which have to vanish. To simplify the discussion, we introduce the notation 
Z^'fi where the first superscript i indicates the loop order, and the second superscript j gives the order of the 
pole, i.e., Z« = ELi Z<H)/e>. 
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When removing the pole proportional to 1/e x In 2 we obtain Z^' 2 ^ as a function of Z^ x \ zjj 1 ' 1 ' , Z^ 1 ^ and Z^f'. 

Similarly, the removal of 1/e x In gives Z^^ as a function of Z^' 1 ^ and Z^ '■ Finally, the removal of the pole 

proportional to 1/e 2 x In gives #g as a function of Z^ and . Inserting this into the expression for Z^ ,2 \ 

we obtain 
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Taking into account all these relations, the 1/e™ pole terms are found to be removed with Z v : and Z\, ' which are 
functions of Z^ only, while Zy 3 ' 3 ^ is determined independently of Zj^ . 

Using the known results for Z^' 1 ' and Z^~\ we recover the result ( |l5|) derived before in Ref. [Q. Using the known 



result for Z^' 1 ', Z^ 1 ,, Z^'S*. Z 2 \ we also obtain 
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this coinciding with previous results derived from a renormalization of the two-point functions in Refs. |Tl], |24|] . 
Finally, inserting Z^' 1 ^ into Zy 3 ' 1 ^ and Z v 3,2 \ we have 
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where, as mentioned above, Z v is independent of Z A ' . For e 2 — we recover the three-loop result of the pure </) 4 
theory of Ref. pit]. 



VI. RENORMALIZATION GROUP FUNCTION OF THE VACUUM j v 



The critical behavior of the renormalization constant Z v is characterized by the finite renormalization group function 
7„ defined by the logarithmic derivative [^l| 
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Using standard methods p5[ , 7^ can be extracted from the simple pole terms of Z v , whose residue will be denoted 
by Zl 1] , as pi 
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The results of the last section for the Zv yield 
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VII. VACUUM ENERGY DENSITY 



In the previous section, we have focused on the removal of divergences, thus fixing the renormalization constants. 
Since the three-loop integrals /3a — I$j are known to zeroth order in e°, we can also determine the finite vacuum 
energy density of the symmetric phase of the Ginzburg-Landau model. Up to a negative sign, it is given by the sum 
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of the vacuum diagrams. Having in mind the application of our result to phase transitions in three dimensions, we 
give here its e expansion. We must calculate the one-loop diagram up to the order e 2 and the two-loop diagrams up 
to the order e. 

The general form of the L-loop result has an e expansion of the form 
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where we have assumed that e 2 and g are of order e, which is correct at the fixed point relevant for the neighborhood 
of the phase transition. The expansion coefficients are 
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VIII. CONCLUSION 



With the help of dimensional regularization and the modified minimal subtraction scheme MS we have computed 
the vacuum energy density in an e expansion up to three loops for the symmetric phase of the Ginzburg-Landau 
model. Further, we have determined the renormalization group function of the vacuum j v , which is the same above 
and below T c . Both quantities will be needed for the calculation of the amplitude ratios of the specific heat at the 
phase transition of the model. 

To arrive at the final goal of deriving universal amplitude ratios for the specific heat above and below the phase 
transition, we must perform a similar calculation also in the ordered phase below T c . Such a calculation will be 
complicated by a proliferation of Feynman diagrams, the appearance of more mass scales and infrared divergences for 
N > 2. Fortunately, the amplitude ratio of the specific heat does not require knowledge of the full effective potential, 
but only its value at the minimum, whose evaluation is simpler and will be given in future work. 

Certainly it is hoped that a higher loop effective potential describing the phase transition will give us specific 
information on the nature of the superconductive phase transition, in particular on the value of the Ginzburg parameter 
at which the transition becomes tricritical. 
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L,m,n2,n 3 



1,0,0,0 
2,0,0,2 

2,1,0,0 

3,0,0,4 
3,0,1,2 

3,0,2,0 
3,1,0,2 

3,2,0,0 




TABLE I: Relevant one-particle irreducible vacuum diagrams W^'™ 1 '™ 2 '™ 3 ^ and their weights through three-loop order ol 
the O(N) Ginzburg-Landau model, where L denotes the loop order and ni,n2,ri3 count the number of g, e 2 and e vertices, 
respectively. 



